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ABSTRACT 


Many applications in optics, such as the diffraction theory 
of optical aberrations, involves the analysis of a function 
defined on a rotationally symmetric system, with either a circular 
or annular pupil. In order to numerically analyze such systems 
it is typical to expand the given function in terms of a class of 
orthogonal polynomials. Because of their particular properties, 
the Zernike polynomials are especially suited for numerical 
calculations. We develop a recursive algorithm that can be 
used to generate the Zernike polynomials up to a given order. 

The algorithm is recursively defined over J where R(J,N) is the 
Zernike polynomial of degree N obtained by orthogonal izing 
the sequence r*^, r^ + 2, . . . , r^ + 2N over (e, 1) . The terms 
in the preceding row - the (J-l) row - up to the N+l term is 
needed for generating the (J,N) th term. Thus, the algorith generates 
an upper left-triangular table. This algorithm has been placed 
in the computer with the necessary support program also included. 
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INTRODUCTION 


An arbitrary function W(r / 0) / such as a wavefront error 
function over a circular or annular region, can be expanded in 
terms of an orthonormal series of orthogonal polynomials. If 
W is defined over a circular or annular region, it is convenient 
to expand W in terms of the Zernike polynomials, Z n ^. This 


W = 


n,l * - 


. A nlV. 


It can be shown {3} that Z ^ = R l(r)e^® where R„^ depends 

n n n 

only on the radial coordinate and e*^® depends only on the 
angular coordinate. Also, 1 is the minimum exponent of the 
polynomials Z n ^ and R n ^ and the numbers n and 1 are either both 
even or both odd. The radial polynomials R n "^ are of degree n and 
satisfy the relation 



If we write, using only the real part, 


W (r, 0) 





Sin j 0 
Cos j 0 


the following properties are satisfied. 

1. The Zernike polynomials are invar ient in form with 
respect to rotations of axes about the origin {3}. 

2. The Zernike polynomials are easily related to the 
classical aberrations {4}. 
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3. The function W(r,0) is usually found as a best least- 
squares fit to a collection of data points. Since the Zernike 
polynomials are orthogonal over an annular region, the well-known 
minimum-error property of Fourier expressions shows that each 
term 

Sin J 0 

A nj Rj jn+2 Cos J 0 

also represents individually a best least-squares fit to the 
data. Thus, the average amount of each term is given by the 
magnitude of that term, without the need to do a new least- 
squares fit. 

Because the Zernike polynomials are being applied to an 
increasing number of physical problems {1,2, 5, 6}, there is an 
expressed interest in being able to generate the Zernike polynomials 
up to a given order. In this paper, we develop a numerical method 
due to Tatian {8} for r ancsrating the polynomials R^n+J over 
an annular region. A discussion of the derivation of the 
algorithm is presented in section 2 and a general discussion of 
thv computer program that was written to facilitate this 
algorithm is diswussed in section 3. 

Bhatia and Wolf {3} have shown that the polynomials 
R J 2n+ j (*) are obtained by orthogonal! zing the functions 

J J+2n 

r , r , . . . , r over the interval {e,l}. The constant 

e represents the inner radius of the annular region and 1 
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represents the outer radiut. Thus e = o represents a circular 

region. We convert this into the associated problem of 

orthogonalizing 1, u, . . . , u n over £e 2 ,lj with the weight 
i 2 

function u J by substituting u = r . This follows from 

A 


<■ 


J+21 J+2k 


s > 




r J + 21 r J + 2k rdr 




r 2J + 21 + 2K rdr 


(r 2 } J+l+k rdr 



u 


j+l+k du = i 


r 


J 


k . j 


u u (u J du) 


First, let us consider the case J = c. Then 
R® (u) is obtained by orthogonalizing the sequence l,u, . . . , 
u n over fe.lj with the weight function 1. This however is a 
Jacobic problem over the shifted interval fe 2 ,ll. Hence, 

» a J 


(1) R 2n (r) - n 


R- (u) * pi°' 0) (2) 


r 2 -V - 1 

TT7 


where p ^°'°^ is the Legendre polynomial of degree n. 
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Now, let us consider the case when J ^ o. We will show that 
R J 2 n+ 2 can b e obtained by a recursive method over the variable J. 

Suppose we have solved the problem of obtaining Q ^ ^ (u) i.e. 

n 2 J— 1 

orthogonalize l,u, . . . , u oyer {e ,1} v/ith weight function u 

and we want to obtain Q J n (a); orthogonalize l,u, *” u° over 

[ 2 - j 

e , lj with weight function u . We obtain a relationship between 
these 2 polynomials by making use of a special case of the following 
theorems { 7 } . 

Theorem 1. Let {p^ (x) } be the orthonormal polynomials associated 
with the distribution da(x) on the interval (a,b). Also, let 


p(x) = c * (x - x 1 ) * (x - x 2 ) . . . (x - x 1 ) 


be nor. negative in this interval. Then the orthogonal polynomials 

{q (x) } associated with the distribution p(x) da(x) can be 
n 

represented in terms of the polynomials p n (x) as follows: 

I Pn (x) Pn+l (x) *Pn+l (x) I 


P (x ) q n (x) 


Pn(x.j) Pn+1 (x^) ’**Pn+l(x^) 


Pn(xl) 


Pn+1 (xl) 


*Pn+l (x ) 
e 
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Theorem 2. The following relation holds for any 3 consecutive 
orthogonal polynomials: 

(y \ 

P„ <x) - (A „ X+ V Pn-X <x) " C n ’ P„-2 


of the highest coefficient of p n (x) is denoted by k n , we have 


n 


= ^n 


Jc . 
n-1 


and C 


n 


n 


n-1 


Theorem 3. 

P Q (X) P Q (y) + P 1 (X) P 1 + ‘ * * p n (x) p Jy> 


n * n 


= k 


n+1 


P n+1 (x) p n ty) . p n (x) p n+1 , y) 


x-y 


using Theorem 1 with p (u) = Q J ”^ (u) , q (u) = (u) and p(u) 

J n n n n 

we have 


2.1 u rt °(u) =J-1 


(o) 


n 




0n J ' 1(o) 


Cl <»> 


Cl 1 <°> 
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2 

We note that this formula becomes undeterminate for u = r = o. 
But we need the values of at u = o to calculate the next set 
of polynomials, so expand equation 2.1 and use Theorem 3 to get 


.J-l 


(2) £T (u) = n 


n 


(1-e 2 ) 


J-l 

V o) 



Q 


J-l 


(o) 


flj-hu, 


1=0 


J-l 


h. 

i 


where 


(3) h^ = \ 
n 2 


(Q^ (u) } 2 u J du 


is the normalization constant for (u) . By substituting the 
expression in equation 1 for (u) into equation 3 and using 


Theorem 3 , we obtain 


(4) h. = -1 


1-e 


7- 


h ^ 1 

— 1 


O'T“ 1 (0) 


converting back to the variable r, we have 
(5) R 2n+J <r) * rJ Q n J (r2) - 
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3. Some comments about the computer program. 


The outline for this computer program is; 
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For a given value of N, the following table shows the values 

of and H J that are generated, 
n n 


TABLE 1 



NOTE: Since is dependent upon knowuing f each of the 

above rows are shorter by 1 entry. 
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The polynomials Q 
namely. 


o 

n 


are directly computed from forumla 


( 1 ) 


Q° (u) = p (o,o) (au+v) 
n n 

where p n ^° r °^ is the Legendre polynomial of degree n, a=2/(l-e^) 

2 2 

and v = ~(l + e)/(l**e). Likewise, the constants 

H° are directly computed from (3); namely. 



NOTE; The above integral is computed in closed form. This is 

possible because Q° (u) is a polynomial of degree n. This is 

done via a call to SQPOLY (square the polynomial) and a call to 

INTGRL (find the integral of a polynomial) . 

Once the first row is knovn (Q 1 ^ and H J for J = o) , the 

n n 

recursive algorithm can then be used to compute each succeeding 

row (Q^ and for J = J q ) . The results are then printed via 

J 2 J 2 

a call to RJN (compute R 2n+J = r (r ) • 
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4 . Example 

A computer run for n = 3 is shown below. 


ENTER N- LARGEST R(0,N) DESIRED 
?3 

R(0,0) 

I.ww 

**************** 

R (0, 2 ) 

NORM: 3.83767639 

-.5487 + 1.0000*R** 2 

**************** 


R(0, 4) 

NORM: 16.46613630 

.2331 - 1.0973*R**2 + 1.0000*R** 4 

************************ 


R (0, 6) 

NORM: 71.94691818 

-.0981 = . 7809*R** 2 - 1.6460*R** 4 + 

************************* 


R(l» 1) 

NORM: 1.90923153 

1. 0000*R** 1 + 

**************** 


R(l, 3) 

NORM: 8.32586779 

-6724*R** 1 + 1. 0000*R**3 

********** 1 \ ****** 


R(l, 5) 

NORM: 35,89625267 

,3191*R** 1 - 1 . 2252*R** 3 + 1.0000*R** 5 

************************* 

R (2 2) 

NORM: 2.32829044 

1. 0000*R** 2 + 

************** 


1. 0000*R** 6 
******* 
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R(2, 4) 

NORM: 11.84460583 

- . 7897*R** 2 + 1.0000*R** 4 

****************** 


R(3,3) 

NORM: 2.6873584 

1.0000*R** 3 + 

*************** 


* STOP * 0 

A COPY OF THE ABOVE PROGRAM CAN BE OBTAINED FROM THE AUTHOR. 


Xlll-11 



PEFERENCES 


ORIGINAL PAGE IS 
OF POOR QUALITY 


1. Baraket, Richard, Optimum Balanced Wavefront Aberrations for 
Radially Symmetric Amplitude Disturbances; Generalizations 
of Zernike Polynomials , J. Opt. Soc. Am. 70', (1980). 

2 . Bezdid'ko, S. N., The Use of Zernike Polynomials in Optics, 
Sov. J. Opt. Techno 1. , 41, No^ 9, Sept. , (19*74) . 

3. Bhat.ia, A. B. and Wolf, E., On the Circle Polynomials of 
Zernike and Related OrthogonaT Sets, Proc. Camb. Phil. Soc., 

tcttitssitt 

4. Nigboer, B. R. A., The Diffraction Theory of Aberrations, 
Physica (Utrecht) 10, (1947) . 

5. Noll, R. J., Zernike Polynomials and Atmospheric Turbulence, 

J. Opt. Soc. Am. 88, (1976) . 

6. Sumita, H., Orthonormal Expansion of the Aberration Difference 
Function and its Application to Image Evaluation, Japanese J. 

Applied ''PKy?.T - 8\ No. "5," Aug.“Tl9^TT 

7. Szego, G., Orthogonal Polynomials, Am. Soc. Colloquium Pub., 
v. 23, (193977 

8. Tatain, B. , Aberration Balancing in Rotationally Symmetric 

Lenses, J. Opt. Soc. Anu (1974) . 


Xlll- 12 



